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ChaDter 1 
Introduci^ion 


In the modern era, helicopters have established themselves as unique flying vehicles due to 
their flexible manoeuvrability^ and growing utility' in both civil and military applications. 
However, helicopters are basically unstable systems and necessiate the pilots to fly them by 
continuously’ apph'ing the control input. This causes an enormous inci'ease in the pilot’s 
work load and leads to poor riding qualities in response to wind gust and poor visibility. 
In addition, in military’ helicopters this severly aflects the performance of weapon delivery. 
Due to these reasons there is a strong demand for designing an Automatic Flight Control 
System(AFCS) for helicopters to improve their stability and control characteristics, and 
to ensure safety' in flight as well as good flying performance. In short, the objective of 
any AFCS is to provide the pilot with a means to control the vehicle safely and effectively 
throughout the flight envelope, ie, to provide good Handling Qualities(HQ). The design 
and development of AFCS has evolved over the peist four decades from simple mechanical 
to complex sophisticated digital fly’-by’-wire flight control systems. All rotorcraft today 
employ’ some form of feedback control. This may be in the form of an autopilot, a limited 
authority control and stability augmentation system, a full authority fly-by-wirc system 
or even a fully integrated centralised vehicle management system. Typical functions 
provided by’ such an AFCS are 

• Control and stability augmentation for superior Handling Qualities. 

• Automatic flight trajectory and forward speed control. 

• • Wind gust response reduction for extended service of life; ride quality enhancement 
for increased crew and pilot comfort. 
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• Aerodynamic and propulsion performance optimization for reduced operating costs. 

• Flutter suppression and failure reconfiguration for increased safety. 

Since it is very difficult to incorporate all the above functions in an AFCS reliably, most 
of the rotorcraft are fitted with AFCS which can perform the basic functions such as ( 1 ) 
Stability and Control augmentation. (2)Autopilot functions such as regulation of attitude, 
speed and flight path. (3)Signalling or annunciation of mission cues. The most important 
requirement is that the AFCS must provide the desired functions at all weather conditions 
without compromising safety of flight. 

Over the last two decades, .AFCS requirements have been in a state of transition. 
Rotarywdng transmission SN'stems have become more demanding, requiring vehicle man- 
agement systems capable of conducting highly aggressive mission under night or adverse 
weather conditions. Also the performance needs of the modern helicopter requires excep- 
tional agility, manoeuvrability and handling characteristics. This promotes the need for 
control systems of increasing complexity in order to provide acceptable HQ throughout 
the operational environment. The implementation of highly augmented AFCS, new con- 
troller types and pilot information systems will give the designer more flexibility to alter 
the response characteristics of the overall helicopter system and to tailor the desired flying 
qualities. 

On the theoretical front, powerful new techniques for design of multivariable control 
systems such as factorization theory, H^o control, ji synthesis or LQG based Hoo control 
have emerged in the last ten years. These have been succesfully applied to a variety of 
control problems in engineering. In particular,, the popularity of the Hrx, control meth- 
ods [1, 6, 7] are primarily attributed to features such as (i) possibility of combining 
multivariable frequency response and sensitivity properties; (ii) guidelines for transform- 
ing realistic design requirements into performance specifications and (iii) availability of 
state space theoretic algorithms and the associated efficient numerical computational al- 
gorithms. 
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1.1 Problem addressed in the thesis 


In this thesis, a limited but realistic and comprehensive design study is carried out, on a. 
helicopter model, using the technicjues of optimization and factorization theory. More 
specifically the thesis addresses the formulation and solution of the following problems for 
a realistic helicopter model. 

, 1. Design of closed loop control laws for stabilit}’’ augmentation satisfying the MIL- 
H-S501A specification, for two distinct flight conditions, namely at hover and at a 
forward speed of 200kph. 

2. Control augmentation compensator design satisfying the MlL-H-8501 A specification 
for both of the above flight conditions. 

3. Robust closed loop stabilization at hover. 

.A,ny rotorcraft is a highly coupled multivariable plant. Hence the helicopter model 
provides an interesting and challenging case study particularly for achieving performance* 
specification of a standard type such as MlL-H-SoOlA. Now performance specification like 
ADS-33 have also begun to be studied from the point of view of control law design. In 
this thesis however, we shall resrict ourselves to the MIL standard. 

1.2* Organization of the thesis 

The chapters in the rest of the thesis are organized as follows. Chapter 2 discusses the 
dynamics of the helicopter, basic flight controls, major coupling effects present in the 
system etc. It also gives an outline of the HQ performance specifications and describes 
how the flight control system design poses itself as an interesting problem for the control 
engineer. 

In chapter 3, a procedure for linear design is drawn up based on HQ requirements. 
For the two flight conditions namely, hover and forward speed of 200kph, two separate 
two degrees of freedom controllers are designed. These include a stabilizing closed loo|i 
controller for achieving closed loop poles as per MIL-H-8501A specification and a forward 
gain controller for augmenting the pilot's input response. Finally, a robust stabilizing 
closed loop control law is also designed for hovev by considering the multiplicative plant 



uncerainty model. This design is carried out using the modern state space algorithmic 
procedure for tlie Hrx, controller design. 

In Chapter 4 the simulation of the final control laws are carried out and the required 
performances are verified, k full discussion of ail related simulations is also carried out. 

Chapter 5 summarizes the conclusions and presents a discussion on the scope for future 
work. 

1.3 Survey of literature 

Various methods are w'ell known for the design of multivariable flight control .systems. 
One of the most popular methods uses classical single input single output(SISO) tech- 
niques to design one control loop at a time. For instance, the Apache AV05 YAH-64 [b], 
was designed using single loop techniques, with proportional-plu.s-integral(PI) configur- 
ation. Eventhough it gave satisfactoty results, this method can be useful for certain 
problems only. Its capability is extremely limited for highly coupled multivariable .sys- 
tems. Furthermore, analysis of multivariable systems with SISO techniciues can often 
give misleading results. There are two main methods of design of control laws for a 
multivariable system namely, Riccati-based methods such as and structural sin- 

gular value. Eventhough the choice of the design method is not clear, the Riccati-based 
solution methods have a broader acceptability. Doyle et al [3], has pioneered the state 
space approach to H 2 and control problems which has found many applications in 
recent times. The reference [15] discusses the control law design using Hoc methods, of 
a multivariable controller to give decoupled pitch and depth control of a submarine. In 
addition to the basic multivariable design, the paper also considers the development of a 
controller conditioning strategy to ensure that performance is maintained in the face of 
actuator saturation. The simulation of controller and desaturation strategy in thi.s paper, 
using the Hoc techniques shows very satisfactoiy performance. 

The application of multivariable design techniques such as Direct Nyquist Array and 
Characteristic Loci method, to achieve decoupled control of the longitudinal dynamics of 
a Vertical/Short Take-off and Landing{VSTOL) aircraft is described in reference [15]. II. 
was demonstrated that the characteristic loci method was difficult to apply succesfully 
in this application. Both performance and robustness were not satisfactory. Work is in 
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progress to apply methods to the VSTOL problem. 

In reference [15], a design example of how //^, optimal control can be scheduled using 
a controller switching strategy for a VSTOL aircraft is given. The linear design makes 
use of an open loop singular value, shaping method. The technique used assigns each 
linear controller design to a portion of the operational envelope. 

Inorder to investigate the handling, control and display requirements for future VSTOL 
aircraft, the UK Ministry of Defence is sponsoring a research programme called Vectored 
thrust Aircraft Advanced flight ControI(VAAC). The research is extending the database 
used to derive the handling qualities requirements for conventional take-off and landing 
configuration to future STOVL aircraft. The flight control and handling research with 
the VAAC Harrier is given in [12], 

A collaborative research programme has been initiated between Royal Aerospace Es- 
tablishment(RAE) Bedford, Leicester Universit}" and Cambridge Control to build on the 
work done by Yue and Postlethwaite [17]. Yue and Postlethwaite, in [17] presented an 
application of H,x.' optimization theory to the design of helicopter control system. The con- 
troller designed was evaluated at RAE, Bedford and the results showed that this controller 
conforms to the military handling qualities specification. According to the literature, the 
piloted simulation results were excellent showing that it was feasible to use an con- 
troller for improving the handling qualities of a typical high performance helicopter. 

Reference [11] presents the synthesis and analysis of multivariable robust control sys- 
tems for a UH-60 Black Hawk helicopter in forward flight based on the ADS-33 specific- 
ation and the structured singular value analysis. Several robust controllers with different 
configurations were constructed to determine which configuration provides the best overall 
HQ performance. 

Despite the emergence of the extensive control design research programmes for heli- 
copters in the Western countries, similar published work in India is almost absent. This 
thesis is one of the modest attempts at beginning an indigeneous effort in control system 
design program for helicopters. 



CIiaDter 2 


ne’icopter Dynamics 

2 . 1 Intro ducti on 

The design of control system requires a thorougli appreciation of the dynamics of the plant 
to be controlled.. Flight control system is no exception. So, for designing an AFCS for a 
helicopter and evaluating its performance, it is essential to understand the fundamentals 
of helicopter dynamics and its control. The basic function of AFCS design is to provide 
a stable and well behaved platform for undertaking the required mission. Helicopters use 
rotating wings (rotors) to produce lift, propulsive power and control. The consequent 
versatility of flight envelope from vertical take- off/hovering to high speed flight explains 
their wide spread use and continuing efforts to improve their performance. 

2.2 Basic flight controls 

The flight control of a helicopter is achieved by providing the following four control inputs. 
They are (i) Collective pitch of the main rotor, denoted as Oq, for increasing thrust or 
lift of the main rotor, (ii) Longitudinal cyclic pitch, denoted as Oi^, for tilting the thrust 
vector in longitudinal direction. (iii)Lateral cyclic pitch, denoted as Ou, for tilting the 
thrust vector in lateral direction. (iv)Tail rotor collective pitch, denoted as 6ot, for yaw 
control. See Fig.(2.1). These inputs are used as control inputs in any fly-by-wire system. 
Their functions are described in more detail below [4]. 

1. Main rotor collective By collective input, the pilot varies the average pitch angle 
of all the main rotor blades and thereby provides a method of controlling the total 
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Figure 2.1: Control System of a Conventional Helicopter 

rotor thrust and lift. The rotor speed is controlled by an independent engine control 
system, which governs the rotor speed. The engine contol system should act in co- 
ordination with main rotor collective because any increase in Angle of Attack(AOA) 
will also increase the drag on the blades tending to slow down the rotor. Since tlic 
rotor speed should remain constant, there should be a proportionate change in power 
to compensate for this change in drag. This control is, carried out by the engine 
control system. 

2. Cyclic piich controls Cyclic pitch controls cause the rotor disc to tilt and change 
the direction of lift force without changing its magnitude. The rotor is tilted by 


( 









changing the pitch of the rotor blades individually once during the cycle of rotation. 

3. Tail rotor collective The main iMirpose of the tail rotor is to counteract the 
main rotor torque. In the absence of the counteracting torque, the helicopter has 
completely unstable attitude. This control is used differently between the hover and 
forward flight. In hovering it is used to rotate the helicopter alsout a vertical axis 
and controKheading. In forward flight it is used to ensure co-ordinated flight, ie. 
to prevent side slip. 

Mathematicalh% the control input to k-th blade of the main rotor can be expressed as 

Ok = 0q + Oic cosii't; + Ou sm4'k 
where = ^(A: — 1), k = 1 ... N 

4'k represents the azimuth location of the k-th blade and N is the total number of blades 
in the main rotor system. 

2.3 Major coupling effects. 

In almost all flight conditions, the uncompensated helicopter is unstable so the pilot is 
required to make continual small adjustments to maintain stabilit}^ even in a trimmed 
state. In addition a number of interactions are present. 

1. Main rotor collective - Tail rotor A change in main rotor collective always requires 
a change in tail rotor collective to balance the change in torque and a corresponding 
change in power is recjuired to maintain rotor speed. 

2. Cyclic -Collective When longitudinal or lateral cyclic is applied to tilt the rotor 
disk, to increase forward speed or to roll, to initiate a turn, the pilot must increase 
the collective to maintain height. This is because tilting of the thrust vector reduces 
its vertical component. 

3. Longitudinal Cyclic - Lateral Cyclic Due to gyroscopic effects of the rotor system, 
pitching and rolling degrees of freedom of the helicopter are coupled. The moment 
of inertia of a helicopter in roll is about one fourth of that in pitch. Consequently, 

. a longitudinal cyclic input requires significant correction with the lateral cyclic. 
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2.4 Equations of motion 


The equations of motion of a rigid helicopter can be written in the body fixed coordinate 
system. The rigid body has both translational and rotational motion. 'I’he dynamical 

t 

equations of motion of a helicopter is given by [10] 

lih — -{iv-hq - Vhr) + XjMa - g sinO 
Vh = —{uhr — lOhp) + YfMa + g cosO siiicf) 
u'h = —{vhP — Uhq) + ZjMa + g cosO coscp 
IxrP — (/jfjr ~ + /x-(?’ + pq) + L 

hyq = ihz - Ixr)rp + Ixz{r^ -p^) + M 
hzr = Uix- Iyy)pq+ Ixz{p-qr)Y N 

( 2 . 1 ) 

0 = p + q sin 0 lanO + t coscb ianO 
9 — q cos 4 ) — r sind) 

4 ' = q sincj) sec9 + cos(j) seed 

The states of the sj'stem and control inputs chosen are 

= [u/,, Vh, lOh, p, q, r, 6 , 4 >, xf}] 

9i3, Oic, ^lo(] 

where Uh, Vh, n'h are the translational velocities along the three orthogonal directions of 
the fuselage fixed axis system, p, q, r are the angular velocities about the same axis and 
6 , 4>, t/) are the Euler angles, defining the orientation of the body axes relative to the 
earth. The control vector has four components: main rotor collective(tlo)i longitudinal 
cyclic(0is), lateral cyclic(0ic) and tail rotor collective(^o<)- Now, Ecin.(2.r) can be written 
in the form 

i = F{x.u) (2.2) 

where F is a nonlinear function of x and u. The solution of Eqn.(2.1) depends upon 
the initial conditions of the motion state vector and the time variation of the vector 
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function F(x,u), that includes aerodynamic loads, gravitational forces, inertial forces aiul 
moments. The trajectory can be computed using any of a number of different numerical 
integration schemes, but gives little insight in to the physics of flight behaviour. We need 
to turn to analytic solutions to obtain a deeper understanding between cause and effect. 
But the scope for deriving analytic solutions of general nonlinear differential equation is 
extremelj" limited. So we go for the linearised approximation of Eqn.(2.2). The essence 
of linearisation is the a.ssumption that the motion can be considered as a perturbation 
about a trim condition. The function can usually be expanded in terms of the motion 
and control variables. 

Similarly, the response can be written in the forfti 


X = x^ + Sx (2.3) 

where .t^ is equilibrium value of the state vector which is referred to as trim and 5x is 
the perturbation away from trim. By doing the linearisation of Eqn.(2.2) at Xg and trim 
control Ue = 0 the following equations are obtained where x denotes the perturbation Sx 
at Xg. 


X = ,4.r + Bug 


(2.4) 


where the state matrix A is given by 

A = 


and the control matrix B is given by 

B = 


C^\ 

\dx ) 




Thus by performing the linearisation of Eqn.(2.1), we will get the system and control 
matrices as given below. Denoting ^ = AT and with derivatives written in semi- 
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normalised form, ie, Xu = XujMa- 


A = 


■Vu 
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0 
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0 

4 - kjQt 

0 

, 1 , ^ 

0 

sinO^ lan0^ 
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QoScc©, 

— 

0 
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-J 

« 

0 

cosd>^ tan©^ 


K 

N'q-k^R,-k^P, 

0 

K 

N'p-kyQ, 

0 
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"^00 ! 

^0„ 



o • 

( 


i Ze,j 

1 ^or 

^^00 ! 


i "0. 

: ^eor 

‘O"!' 

0 r 0 

r 0 


>'0u L 

i. ^'eor 

; "^0 i 


'L'^o; : 

0 1 

0 

! 0 

; 0 


^6u ! ^0. 
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The coefficients of the A matrix are called stability derivatives while that of the B matrix 
are called the control derivatives at trim. Eqn.(2.4) is the fundamental linearised form 
for describing the stability and response of small motion about a trim condition. 

2.5 Handling qualities 

Handling qualities denote the performance of the helicopter (ie, vehicle’s response) to 
pilot’s commands. The original definition of handling qualities by Cooper and Harper 
is [2]: 
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Those qualities or characteristics of an aircraft that governs the ease and precision with 
which a pilot is able to perform the tasks required in support of an aircraft role. 

Goodness, or qiiality, according to Cooper-Harper, can be measured on a scale span- 
ning three levels as shown in Fig.(2^). Aircrafts are required to be Level 1 throughout 
the Operational Flight Envelope(OFE). Level 2 is acceptable in failed and emergency 
situations but Level 3 is considered as unacceptable. 


adequacy hr 
selected task or 
required operation 


aircraft 
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task Pitot 
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to mintmaJ 
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I 
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Figure 2.2: Cooper Harper Handling Qualities Rating Scale 

The requirements of performing a given task drive Handling Qualities R.equirenieiit,s(nQH.s 
and they are given in flying quality standards, like the US Army’s Aeronautical Design 
Standard for handling qualities (ADS-33) and MIL-H-8501A. The criteria in these stand- 
ards quantify responsiveness and sensitivity of the vehicle and lay down mandatory quality 
boundaries on measural)le parameters, which are needed as design targets during design 
and certification stages. • 
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2.5.1 Stability characteristics 

Dynamical stability at trim refers to the abiiitj' of the vehicle to resist itself from diver- 
ging away from trim under small disturbances. The linearized perturbational equations 
are used to obtain the stability characteristics of the base helicopter. The longitudinal 
dynamics of the helicopter has an unstable phugoid (combination of pitch-heave-forward 
translation) which persists through out the speed range. The acceptability of these os- 
cillations, according to the M1L-H-S501A stability requirements in the two flight regimes 
are given in Table 2.1. 


Table 2.1: Summary of MIL-H-S501A Stability Requirements 


Period 

Damping Requirement 


Visual Flight 

Instrument Flight 

< 5 sec 

1/2 amplitude in 2 cycles 

1/2 amplitude in 1 cycle 

5-10 sec 

Atleast lightl}' damped 

1/2 amplitude in 2 cycles 

10-20 sec 

\ Not double in 10 sec 

; Atleast lightly damped 

> 20 sec 

1 

^ No requirement. 

Not double in 20 sec 


2.5.2 Control characteristics 

Control characteristics of helicopter deal with the probleni of evaluating the response of 
the system to a given blade control input. Generally, the study involves analysing the 
response of the system to each of the four control inputs. The main rotor collective is nol 
treated as a short term control. So the response to this input is not considered. 

Pitch, Roll and Yaw response to control inputs 

The total control response is contributed by (i)initial response and (ii)the steady slate 
response. Generally for accurate manoeuvring, it is the initial response that plays a 
major role. The response behaviour is always analysed using the linearised perturbational 
equations developed earlier. 
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Control characteristics or handling quality requirements of a helicopter are specified 
in MIL-H-S501A. These requirements are shown in Fig.{2.3) below. 


Acc*pU6i« UarpirW lV)*oc»pubI« Ymw Axm 



P»lChCootrolPcr»«ref.Tnch Bol Control Y»w Con^ Pow^f/lrch 

Prtch Mon>*n( ot Ro* Monveot ol Y»w Mom#o< oJ 

« f»d s^ r*<Jr»«c^ nci%^ 

n^' k\ in fcn 


Figiiie 2,3: MIL-H-8501A Control Requirements 


These three figures correspond to response in uncoupled pitch, roll and yaw manoeuvres. 
The limits specified are based on pilot's opinion. If the response of the vehicle for a unit 
step pilot input is very slow and small, then the pilot will find the vehicle to be sluggish 
and if the response is too large, then the vehicle is over sensitive. The parameters of the 
two axes correspond to the control power and control damping. These will be explained 
by following the uncoupled motion in any one axis are shown in Fig.(2.3). For example 
the uncoupled roll equation can be written as 


rdL dL , 


(2.5) 


( 2 . 6 ) 




Note : Lp = l^r ; La — 

Lp = Lpj I XX ; Lg^^ = ^&icf 
Its solution is gi\’en by 


p=z 


LeJI. 
Lpf Ixx 


■[1 _ 


(2.7) 
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The handling quality requirements are based on the two parameters representing the 
initial angular acceleration of roll motion for unit control input and the damping in roll, 
viz 

Similar parameters are used for pitch and yaw with the corresponding control moment . 
damping and inertia. 

From the foregoing discussion, it is evident that the control system designer must 
consider the various coupling effects in the system dynamics and also any nonlinearitics 
present in the elements. Due to complexity of the problem, the design of a robust con- 
trol system for a rotorcraft, meeting the prescribed HQ performance specifications under 
anticipated disturbances and uncertainties is a challenging task. 
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Chapter 3 

Design and Development of FLigh" 
Control Laws 

3.1 Introduction 

This chapter is devoted to the development of control laws for stability and control aug- 
mentation for a helicopter. Further the design of a robust stabilizing controller for the 
helicopter model at hover using the multiplicative description of plant uncertainty is also 
discussed. The problem of designing control laws for the AFCS requires modeling of 
l)the vehicle dynamics in different trim conditions and understanding the extent of mod- 
elling uncertainties and 2) a quantitative list of specifications on performance. These were 
discussed in the last chapter. 

3.2 Helicopter model description 

For FCS design, it is necessary to consider the dynamics of other components such as 
the cyclic mixer, the actuator and the mechanical linkage between the pilot sticks and 
the actuator in addition to the vehicle dynamics. At a functional level the closed loop 
system formed Iw these subsystems is described schematically in Fig.(3.i). For each of 
the functional blocks shown above the d^mamic equations are developed below. For the 
vehicle dynamics these are linearised equations at a specified trim condition. 

1. Vehicle dynamics: A detailed description of the vehicle dynamics is given in chapter 
2. For the Stability and Control Augmentation System(SCAS)design, only the cyclic 
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Figure 3.1: Vehicle and Components 
inputs are considered. The equations of motion is given by 

X = .4.r + Bu 

where 


(3.1) 


x^ = [u, c, 10 , p, q, r, 6, 4>] 

t/ = Oic] 

The decoupled longitudinal and lateral dynamics is used in Control Law(CLAW) 
design and is written as 




{:t.2) 


with 


= [“> 0] 

^In — [^Ic] 


for longitudinal axis and 


xit — AitXit + BitUit, 


(3.3) 
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with 


•I’/T = [fN 4 ^], 

uit = [6'ic] 

9 

for lateral axis. The model of the helicopter at hover is given in Appendix D. From 
this data, the matrices in Eqns.(.3.2) and (3.3) for the decopled longitudinal and 
lateral dynamics are as follows. 


Mn = 


Alt = 


-0.0157 

0.0312 

0.00.53 

-0.1703 

0.0048 

-0..3.346 

0.0045 

-0.0172 

0.9390 

-0.0905 

-1.3220 

0.0000 

0.0000 

0.0000 

0.9900 

0.0000 


; -0.2061 - 


Bln 

' -0.0103 ; 

21. .5493 



0.0000 


-0.0001 

-0.0084 

0.0090 

0.1701 

-6.6900 

-5.1620 

-0.0426 

0.0000 - 

5.6710 

-0.56.50 

-0.9319 

0.0000 

0.0000 

1.0000 

0.1012 

0.0000 


(3.4) 


Bit = 


0.1950 

3.0026 

14..3031 

0.0000 


(3.5) 


'Is 


Cyclic mixer: For the cyclic mixer shown in Fig.(3.1) the equations are algebraic, 
ie, 

■ 9u ■ ^ 

The effect of cyclic mixer has already been included in the matrices A and B, given 
for the particular model. 


(3.6) 
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3. The actf;<i|prs: The actuators are modelled as first order lags with time constants 
Ti,r 2 both between 25 and 100ms. Their equations are 

' I _ I 0 ' ^1* : ^ ' ^la?J + ^Un ' 

^ 0\c j i ^ ■“'1/^2 j ■ ^Ic r 0 l/'^2 ■ \ d\cp + ^Ica 

For the given problem , the value of both ti and ra have been taken as 25ms. 

4. Pilot's inputs: The components of inputs Onp and O^p are obtained from the pilot’s 
inputs through the mechanical linkage. The function of these linkages is to scale 
the pilot's inputs by constants. The particular A and B in Eqn.(3.2) includes Uie 
dynamics of the linkage also. 

In order to design the control law, the above models have to be connected together to 
form an appropriate i/p-o/p system. This is done in the later sections. 

3.3 MIL-H-8501 performance specification 

The goal of any rotorcraft control system design is to achieve ” Level 1” handling qualities 
performance. The M1L-H-S501A specifications for handling qualities is listed in chapter 2 
and is reproduced below. As seen in chapter 2, the unstable phugoid persists throughout 
the flight envelope. The acceptability of oscillations in the two flight regimes are given in 
the Table 3.1. 

Table 3.1: Summary of MIL-H-8501 A Stability Requirements 

V, 

, Period Damping Requirements 

Visual Flight Instrument Flight i 

, < 5 sec , 1/2 amplitude in 2 cycles , 1/2 amplitude in 1 cycle 
5-10 sec Atleast lightly damped ‘ 1/2 amplitude in 2 cycles 
. 10-20 sec Not double in 10 sec Atleast lightly damped 

> 20 sec No requirement > Not double in 20 sec ; 

The control characteristics or handling quality requirements of a helicopter, as specified 
in MIL-H-S501A are shown in the Fig.(3.2) ' 
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Figure 3.2: MIL-H-8501A Control Requirements 

3.4 Approach to control law design 

The control law design for stability and control augmentation uses the decoupled longit- 
udinal and lateral dynamics of the helicopter and separate controllers are developed for 
each axis. The closed loop system including the controller should satisfy the perform- 
ance requirements listed in Table 3.1 and Fig.(3.2). But the two types of performance 
requirements ie. the stability requirement and the handling requirement are interdepend- 
ent. Hence the design procedure involves designing the closed loop law satisfying the 
stability requirements first, followed by the development of forward compensator to tailor 
the handling requirements [18]. 

3.4.1 Stabilization 

The design of control laws for stabilization is first considered. For this purpose, augmented 
system model, which is the series combination of the vehicle dynamics and the actuator 
dynamics is required. The state vector, input vector and output vector for the augmented 
model in the longitudinal axis are chosen as, 

xf = [ u. It', q, 0, 0\^ 

9n 



Ilf = 

I/f = [ 9 .^] 

The state vector above consists of vector components xin, of the vehicle dynamics and 
•To, those of the actuator dynamics. The dynamic model, obtained by connecting these 
components, using the data given, is 

X “f" Bauglll ('.’•b) 

y ^'augXi T Dauglll 


with 


-0.0157 

0.0312 

0.004S 

-0.3346 

0.9390 

-0.0905 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 ^ 


0.0000 : 



^avg — 0.0000 

^ 0.0000 
1.0000 


0.0053 

-0.1703 

-8.2440 

0.0045 

-0.0172 

-0.4122 

-1.3220 

0.0000 

861.9700 

0.9900 

0.0000 

0.0000 

0.0000 

0.0000 

-40.0000 



-00 1 00 

; 0 0 0 1 0 


D 


aug 



Stabilization problem: To find a controller with input y and output u such that the 
closed loop system has eigen values in the specified region so as to satisfy the stability 
requirements given in Table .3.1. Using the parametrization formula a family of controllers 
satisfying this criteria can be found(see Appendix A). One such controller is given below. 
The steps leading to this controller are as follows. 


1. Fix a set of 5 distinct complex numbers so as to, satisfy the criteria in Table 3.1. 
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2. Find matrices K and H such that A-BK and A-HC have eigen values at the above 
complex numbers. 

3. The formula for the controller transfer matrix is then given by 

C{s} = X{s)Y-^{s) (3.9) 

where 

^ A-BK 
X{s) = ; 

V 0 ; 

^ A-BK 

y(s) = ; 

^ C - DK I j 

4. The stability margin.closed loop sensitivity function and closed loop eigen values 
are then determined. 

In order to design the longitudinal controller, the complex numbers in step 1 are chosen 
as 

p = -40 -1.4050 -0.3386 -0.2±0.333h‘ 

The matrices K and H are obtained as 

K = ■ -0.0990 -0.0Q58 0.0213 0.0429 0.4713 

’ 0.3992 -0.1703 : 

0.0317 -0.0172 

-H= -0.9337 0.0000 

0.9900 1.4050 - 

’ 0.0000 0.0000 

The longitudinal controller at hover is then obtained as 

f Ac Bc^ 

C{s) = ^ (3.10) 

^ Cc Dc J 
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where ,4c, Bc~ CV. Dc are 


■ -0.3636 

-0.4560 

20.4644 

.34.4671 

170.6.367 ■ 

-0.1194 

-2.1.508 

74. .50 13 

125.8861 

621.7803 : 

^ 0.0000 

0.8774 

-63.29.56 

-105..5043 

-525.4513 ' 

0.0000 

0.0000 

14.8447 

23.0912 

118.2093 ; 

= 0.0000 

0.0000 

0.0000 

0.2468 

0.1039 

-0.1.360 

0.0902 




. -0.0571 

0.9339 




1.8189 

0.7776 




. -0.8317 

-1..5719 




-0.0499 

0.2257 





Cc = , 0.0000 0.0000 0.0000 0.0000 0.3423 : 

Dc = 0.0000 0.0000 ' 

Moving along the same lines, the controller in the lateral direction can be obtained for the 
given model at hover. It is given below. The complex numbers in step 1 are chosen as 

P = . -0.2 ± 0.3.331?' -5.21SS -2.0000 -40.0000 
The lateral controller is 


0.3778 

4.1278 

-181.02.53 

187.2714 

473.1824 

-0.4111 

-4.5483 

66.4548 

-.58.9959 

-149.42.50 

0.0000 

0.4072 

7.4586 

-9.9965 

-23.8.347 

0.0000 

0.0000 

47.8798 

-50.8002 

-127.0600 

: 0.0000 

0.0000 

0.0000 

0.4727 

-2.0315 

^ 1.247 

12.1144 




0.9881 

-3.0068 ' 




-0..3516 

-0.9236 ■ 




-0.8198 

-4.1664 ■ 




' 0.1690 

0.8303 





Cc = 0.0000 0.0000 0.0000 0.0000 1..3898 

Dc = 0.0000 0.0000 - 
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The same procedure is carried out for the helicopter model at 200kph. Tlie model is given 
in Appendix D. The longitudinal controller has been found out and is given below. The 
pole location chosen in the design criteria is 


p 

= ; -0.2 ±0.33311 

1.4050 - 

0.3386 - 

40.0000 = 

The controller is 

^ -0.2829 

-0.3710 

-26.5326 

14.8712 

446.5495 


: 0.0891 

-0.4112 

5.1943 

-1.6459 

-77.0964 

.4, = 

0.0000 

0.6779 

-27.7090 

12.8070 

432.6279 


; 0.0000 

0.0000 

28.0092 

-14.7179 

-458.3225 


; 0.0000 

0.0000 

0.0000 

0.1179 

-1.5189 


^ 2.0253 

3.7937 





' 0.9013 

0.5828 




Be = 

0.2947 

0.8890 1 





-2.8761 

-3.1318 ' 





0.0984 

0.1166 ^ 




Cc = 

0.0000 

0.0000 0.0000 0.0000 

1.7695 ^ 


De = 

0.0000 

L 

0.0000 ^ 
j 





Lateral dynamics of the model at 200kph is stable. Its eigen, values are 

-0.6894 ± 2.3547f -5.8189 -0.0416 -40.0000 
So there is no need of designing the control law in the lateral direction. 

Now the closed loop stabilizing control law for the whole coupled system at hover has 
been found out using the same design strategy. The controller has been found out and Is 
given below. The pole location chosen in step 1 is 
P=,-40.0 -40.0 0.5 -2.0 -0.5 ±1.271 -0.3 ±0.41 -0.2 ± 0.351 i 
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The controller is Ac= 


: -0.27 

1.07 

-19.30 

-48.22 

195.90 

169.25 

-84.25 

-67.44 

-503.75 

587.69 

1 

p 

o 

-0.47 

17.44 

-84.25 

328.24 

365.19 

—35.13 

-92.78 

380.21 

-184. 7'. 

0.01 

0.07 

o.-^g 

-8.18 

17.41 

17.30 

1.54 

-5,44 

1.05 

2,15 

0.01 

O.OS 

0.72 

-6.65 

10.26 

9.28 

3.51 

-2.07 

2.17 

-0.97 

0.00 

0.00 

-0..52 

4.43 

-22.49 

-24.78 

5.26 

7..32 

-3.42 

-11.02 

! 0.00 

0.00 

0.68 

3.39 

-11.17 

-18.01 

1.07 

-6.93 

2.49 

-8.35 

■ 0.00 

0.00 

0.00 

0.00 

5.57 

4.48 

-2.97 

-5.61 

-1.10 

1.56 

0.00 

0.00 

0.00 

0.00 

3.42 

-7.30 

-4.01 

-37.14 

0.74 

2.64 

■ 0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

1.64 

-2.4 

-0.81 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.12 

0.16 

2.39 

-0.42 



; 44.71 

-7.22 

2.-32 

-88.26 







■ —34.81 

3.15 

-2.46 

15.51 







0.30 

0.31 

0.86 

-0.93 







- 0.32 

2.07 

0.89 

-1.00 






Be = 

0.27 

0.59 

0.53 

0.76 







-0..54 

-16.84 

-0.30 

-1.51 







-0.08 

-6.09 

-0.35 

-0.70 







-1.29 

-53.57 

0.10 

0.49 







0.01 

2.25 

-0.12 

-0.58 







0.00 

0.24 

-0.01 

0.03 






Ce = 

0 0 0 

0 0 0 

0 0 

1.90 0..30 

De 

" 0 0 

0 0 




'000 

0 0 0 

0 0 

-0.10 8.30 

0 0 

0 0 



3.4.2 Handling qualities criterion 

The next step is to design a forward compensator to achieve the desired response qual- 
ities. For this the closed loop equations have to be formed. The transfer function of the 
closed loop system with 6ic as the input and y as output is given by 


T 

1 yp 


^ Act Bel \ 
V Cci 0 j 


(3.11) 


2.5 



wlicre 


t 


.4.-; = i 


|~ 1 nng ^ aiig 


L 


^ ting 


^aug^’c 

Ac 


Bel = 
Cel = 


L 


I! 115 

0 


c. 


'“.7 
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I'n orcier to satisf}’ ihe specifications on handling qualities, a forward coiniKMisai or 
designed so as to provide additional i/p summed up at the acuator i/p. The gaiii of 
this forward compensator is then chosen to satisfy the HQRs of Fig. 3. 2. d'hc rclevaiii 
equations are 


1. (flosed loop equations : 


Xci — Acl-l'cl + BclVp 
■ y = CrlXel 

2 . Media 11 leal mixer : 

= A'p?/ 


where y is the pilot input. 

3. Foinvard compensation block : 


.r = F;/ 


4. Summation of inputs : 


Vp — {hp + F)i] 


5. Compensated system : 


.I'd = AciXd + Bciihp + A )// 

y = Cct-i'ci 
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The above formula of tlie dynamical system from pilot input r; to the output y then involves 
the gain block F which can be adjusted, thus meeting the response requirements without 
affecting the stability requirements met b}' the closed loop controller. The uncoupled 
equation for pitch in the closed loop system formed, for the problem at hover can be 
written as 

<7 = — i.322<7 + 0tf 

Now, a scalar /i is chosen as fi = 0.2 such that 

q = -1.-322^ + /iu (3.12) 

belongs to the acceptable region of Fig.(3.2). Similarly for the model at 200kph, the 
equation for pitch is 

q = —1.7540(7 + Ou 

for which the scalar fi chosen is fi = 0.2 inorder for the control reejuirements to be met. 
Similarly the uncoupled equation for p and r at hover can be written as 

p = — 5.1620p + 0if 
r = — 0.93r9r + Ou 

As in the longitudinal case, we have to find out scalars /2 and fz such that 

p = — 5.1620p + / 2 U 
r = —0.9319?’ + fzu 

meet the control specifications. For this fz is chosen as /a = 2, but fz cannot be found 
out, since the scale given in Fig.(3.2) for p is not sufficient. For the model at 200kph the 
equations for p and r are 

p = —5.592/) + Ou 

r = — 1.4167’ + 0u 

Here also the scalar fz can not be found out since the scale given is not sufficient, fz 
has been found out to be fz = 2, for the control specifications to be met for the model at 
200kph. 
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3.5 Design of control laws using PI techniques 


Let us try to analyse the problem of closed loop stabilization using a PI controller. The 
output signal produced by the i’l controller consists of two terms, one is proportional 
to the input signal and the otlicr proportional to its integral. For longitudinal hover tin* 
output signal of the controller is 

C{s) = A'l / ' ^ : dt + A 2 : 

0 e 

A'l and A '2 can be calculated so that the resulting closed loop system meets the MIL-H- 
8501 stability reciuirements. The values of A'l and A 2 that makes the closed loop system 
stable are given in the Table 3.3. 


Table 3.2: Values of Ah and Ah for the PI Controller 

' Ah Ah 
0.11 0.33 ’ 

0.33 LOO 
0.56 . 0.56 
0.78 0.1 
1.00 0.78 ' 

, 1.00 LOO 


But in PI methods the compensator for arbitrary pole locations may not exist, unlike the 
case for state feedback compensator. Also thePI compensator may not give satisfactory 
response for the whole coupled system. 

i 

3.6 Robustness 

No nominal plant model can emulate a physral plant perfectly. Thus in the design of 
control system, we are faced with the uncertinty of the model of the plant, originated 
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from various sources such as the neglected nonlinearity of the plant dynamics, the sim- 
plification of the plant model for the purpose of controller design, the variation of plant 
parameters during the operation and so on. For instance, the major uncertainties as- 
sociated with a rotorcraft model are due to, neglected rotor dynamics and unmodelled 
rate limit nonlinearities. Also, the aerodynamic coefficient vary with flight environment. 
However, the controller (designed for the simplified linear model) should maintain closed 
loop stabilit}^ for the actual system in the face of such perturbations. This is the main 
concern of the robust stability problem. As shall be seen in the subsequent sections, this 
problem is elegantly addressed in the /foo framework. 


3.6.1 Robust stabilization as a standard Hqq optimization prob- 
lem 


In this section we pose the problem of robust stabilization as a standard 4-block problem. 
It is well known that, by making use of the stable coprime factorization theory [14], the 
problem of robust stabilization can be formulated as an H<x, model matching problem. 
In this problem a stable transfer matrix R is determined such that the error j] Ti — 
T^RTz l|oo < 1 where Ti, T 2 , Tz are transfer matrices obtained from the plant and the 
uncertainty profile. However an algorithmic solution of the model matching problem has 
many drawbacks, particularly for the multivariable case. For this reason, it is desirable 
if a model matching problem can be formulated as a standard 4-block problem. This is 
carried out below. 

Robust stabilization at hover 

First we would like to determine whether a robust stabilizing controller can be developed 
at hover such that this controller also stabilizes the model at full speed. This analysis is 
carried out using the multiplicative description of uncertainty. Let Po denotes the transfer 
matrix of the dynamics that corresponds to hover. It is given by 


PoiM 


P.W^) A 5 '' 
F,(ju.) . \C D I 


(3.13) 


Recall that this is the model obtained in Eqn.(3.4). Let A(s) denotes any stable transfer 
matrix denoting a ball of uncertainty such that [] A ||oo< 1. Then an uncertain set of 
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plants which includes Pq. is represented by the set [9, 16] 

M{PoAV) = {PlP = (/ + AH-')Po. i| AlU<l} (:i.l4) 

1 lii.s iii called mutt iplical ive uiicerlaiiily aixMiiid the iiuniiiial plant model Pq. Here V\’ 
is chustMi Hticli (hid M (’uvt'i's (lit* retniirt'd set of iiiotlelH uf the titieei'tiiiii plant. I'ur the 
present purpose we need to find a transfer function weight W such that if Pq denotes 
the model at hover, then the model at full forward speed is also included in the set M. 
Frequency response design of W is described next. From Eqn.(.S.14), it follows that 

II P - Po It (u.-) ' < II IF Po II (cu) Vcu € [0, oo) (.3.1.5) 


The full order model at 330kph is given in Appendix D. From this data, the decoupled 
longitudinal model P is written as 




PiU<^) 


• W is chosen to be of the simple diagonal form 


^ A b\ 

\C Dj 


(.3.16) 


ai(s) 0 

IF = ^ 

0 (his) ; 

where ai(s) and 02 ( 5 ) are selected such that Eqn.(3.15) is satisfied, ie, 

I Pi{j^ - PiU^') P + I A(iw) - F2(itu) P < I ai{juj)Pi{juj) P + I a2[juj)P2{juj) P 


Using the given data in Eqn.(3.13) and (3.16), ai(s) and 02 ( 5 ) are found out to be 


ai(s) 


02 ( 5 ) 


(10s + l)(0.1s + l) 
4s^ + 2.Ss + 1 
2(0.1s + 1) 

2s + 1 


Robust Stability Condition 

The robust stabilizability condition for the above multiplicative description of the uncer- 
tainty around a nominal plant model Pq is given by the following theorm [14]. 

Theonm- 

Let M{Po,W) and 1F(s) are specified and suppose C stabilizes Pq. Then C stabilizes all 
P in the class M ( Po, IF ) if and only if 


11 [lFPoC(/ + PoC)-*](iu.') Hoc. < 1 Vu; € P (3.17) 
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I he robust stal)ility condition of 17) can be formulated as a 4-l)lock problem as 

given below. I he standard set up is shown in Fig. (3. 3) for which the transfer matrix from 
w to z is [see Appendix B] 

T =n- = Gii — Gi2(/ + G 22 G) ^^21 (3.18) 


\xi 

u 

\ 






G 


z 

>■ 

a 


C < 

I 


Figure 3.3: The Standard Block Diagram 


We can convert ijie robust' stability condition to the set-up of the standard problem by 
<l(dining G so that in Fig. (3. 3). the transfer matri.x from w to z ecpials W PqC[I /u)~'- 
'Fhis is accomplished by taking G'n = 0, GV 2 = — I'FPo, G 21 = /, G 22 = Pq- Hence the' 
■1-block plant corresponding to the robust stability condition is 


with the state space model 


G 


' 0 
/ 



( A Bi B2 ^ 

G = Cl Dll Di2 

u 

^ C2 D21 D22 j 


(3.19) 


While checking the soh'ability of the standard 4-block problem using model matching 
theory[see Appendix C], it was found that || Ti{z} ||= 2.0069, which should be less than 1 
inorder for the robust stability problem to be solvable. Thus it turned out that the robust 


stabilization of the whole flight envelope b}' including the two extreme flight configurations 
as modelled by the above W is not possiljle. Since there exist manyf otlier possible 
functions W which meet the requirement of including P in M in E(in.(3.1d), the above' 
failure to stabilize both hover and full forward s[)eed is not conclusive. However il, ran 
be said that this would involve very comple.x coutrolleM’s and hence' may be impractie’al. 

Norv the robust stal)ilization at hover is carried out with a different VV. ie, 


0 

where erj(.s) and are se-le'cte'd such that l.lu're is nert mudi ele'viatieni be'twe'e'ii 1’ aiiel 

P at low' freejuenries and at high freepiencies there is large eleviation. I’his descriptiern erf 
uncerainty' is infact more realistic since for most dymamic sy'stems the model uncertainty 
begin from the mid freepiency' range while at very low' frequencies the models are good 
approximation to the jriant. For tlie present problem, «i(.s) and rtjl*’) arc chosen as 



TiS + ki 

O.Ls + 1 

T2-S + A’j 
0.1.S + 1 


The solvability of the robust stabilizability condition is checked using the model matching 
theory'. For different tj and T 2 , the value of || Ti(r) || is given in Table 3.2. In order for the 
4-block problem to be solvable, \\ Ti{z) |j< 1. So ti = 0.3, ki = 0.1, rj = 0.3, ^2 = 0.1 
is taken for aj(s) and or2('5)- 


3.6.2 State space procedure 

We now' use the iterative state space algorithm of Glover and Doyle [8], to find the solution 
to the above problem of making || 1| < 1. This requires the condition (1), (2), (3), 

(4) to be fullfillcd. See .^\i)pcndix 13. Hence these are first vertified for the above data. 

For the particidar calculated above and for Pq of Eqn.(3.13), the state space lep- 
resentation corresponding to Eqn.{3.19) is w'ritten as 
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Table 3.3: Variation of || Ti{z) |1 with r and k 


D 

ki 

r-2 

£2 

li 11 

0.02 

• 0.01 

■ 0.02 

0.01 

0.0386 

0.2 

0.1 

0.02 

0.01 

: 0.1290 

0.2 

0.1 

: 0.2 

0.1 

0.3860 

0.3 

0.1 

• 0.2 

0.1 

0.3938 

0.3 

0.1 

; 0.3 

■ 0.1 

0.4563 

0.03 

0.01 

0.3 

0.1 

0.4324 

0.03 

0.01 

0.03 

0.01 

0.0456 

0.3 

' 0.1 

0.03 

' 0.01 

0.1525 


-40.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

-8.24 

-0.02 

0.03 

0.01 

-0.17 

0.00 

0.00 

-0.41 

0.01 ■ 

-0.33 

0.00 

-0.02 

0.00 

0.00 

861.97 

0.94 

-0.09 

-1.32 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.99 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

-10.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

-10.00 

10 0 0 

0 0 

11 

1 0 

0 0 0 

■T 

0 0 



0 0 0 0 0 0 0 
0 0 0 -3 0 29 0 

Ci = 

0 0 0 0 -3 0 29 


0 0 0 1 0 0 0 
0 0 0 0 1 0 0 


Du = Dn = T>2i = . D22- - 

0 0 0 . ■_ 0 1 _ 0 _ 

For the given data rank £>12 = 0. Hence condition (2) is not satisfied. This is due to 


the fact that the control effect is not included in the performance specificatioiffie, robust 
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stability condition). Hence instead of this performance specification, \vc use tlie following 
regularised specification. 


HToC'(/ + PoC)-‘ 
HuC'(/ + PoC')-* 


(3.20) 


clea.rl}" for H i = cl and small e, if a controller satisfies Eqn(3.20), then it also satisfies 
Ec|n.(3.17). The effect of incorporating H i in the performance specification modifies the 
Di 2 matrix. Hence condition {2)is now satisfied. Even when e is not small, the controller 
obtained for the modified specification of Eqn.(3.20) is still a robust stabilizing control- 
ler(which is more conservative). With the performance specification as in E(in.(3.20), tlu' 
4-block plant is now modified as 



Next, condition (3) and (4) are checked for the given data. The first system of condition 
(3) does not have any transmission zeros. For the second system no transmission zero 
is found to be hung on the imaginary axis. Since the transmission zeroes, if they exist 
are at the eigen values of A, a better numerical procedure to check condition(3) and (4) 
is to calculate the condition numbers of these matrices at the eigen values of A matrix. 
Thus it was verified that the given data satisfies both condition (3) and (4). Once the.se 
assumptions on the generalised plant are satisfied, we next check the existence conditions 
of the theorm given in Appendix B. Using the state space algorithm, the suboptirnal 
controller for this problem is given by 


(a b\ 






^vith 


-40.2498 

0.00.55 

0.0385 

42.8005 

7.0336 

3-58.09-52 

-67.8392 

-8.2440 

—0.01577 

0.0312 

-3.20-54 

-3.2860 

30.9651 

29.9435 

-0.4122 

0.0048 

-0.3346 

36.9128 

48.7944 

-356.2358 

-469.3768 

861.9700 

0.9390 

-0.0905 

0.3285 

0.3-571 

-15.8961 

-3.4087 

0.0000 

0.0000 

0.0000 

1.4783 

0.7750 

-4.7164 

-7.4-548 

0.0000 

0.0000 

0.0000 

0.3725 

0.5790 

-3.9325 

—5.-5642 

0.0000 

0.0000 

0.0000 

-0.0195 

1.0569 

0.1889 

-10.-5477 


-0.1360 -0.0960 
0.0074 O.OISI . 

-0.0563 -0.2.554 
B = -0.0061 -0.0045 

-0.0004 -0.0038 . 

-0.0002 -0-0034 
0.0000 - 0.0002 

C = -0.24980.0055 0.0-385 5.6181 -0.0884 0.0198 0.0792 , 

D = 0.0000 0.0000 

Thus the robust stabilizing controller has been formulated as a standard optimization 
problem and a stabilizing controller has been found. 

3.7 Summary of control laws designed 

1. Flight condition : Hover 

• Stability augmentation closed loop law: 

Uc, = Cl ■ 

$ 





for the longitudinal axis with 




Cl = 

ft, 






Iv ft, D„ 



where 







-0.3636 

-0.4560 

20.4644 

.34.4671 

170.6367 


-0.1194 

-2.1.508 

74.5013 

125.8861 

621.7803 


= 0.0000 

0.8774 

-63.2956 

-105.5043 

-.525.4513 


0.0000 

0.0000 

14.8447 

23.0912 

118.2093 


■ 0.0000 

0.0000 

0.0000 

0.2468 

0 . 10.39 


-0.1360 

0.0902 





-0.0571 

0.9339 





= 1.8189 

0.7776 





-0.8317 

-1..5719 





-0.0499 

0.2257 




a. 

= 0.0000 ( 

3.0000 0.0000 0.0000 

0.3423 ' 


Dci 

= 0.0000 ( 

3.0000 


. 


and for the lateral axis 







yc2 = 

r 

C 2 ’ ; 






0 





C2 = 

( A,, Be, 






Be, 

y 


where 







0.3778 

4.1278 

-181.02.53 

187.2714 

473.1824 


-0.4111 

-4. .5483 

66.4548 

-58.9959 

-149.4250 


= 0.0000 

0.4072 

7.4.586 

-9.9965 

-23.8.347 


0.0000 

0.0000 

47.8798 

-50.8002 

-127.0600 


0.0000 

0.0000 

0.0000 

0.4727 

-2.0315 


1.247 12.1144 

0.9881 - 3.0068 

= - 0 . 3.516 - 0.9236 ' 

- 0.8198 - 4.1664 
0.1690 0.8303 

r 

C -2 = 0.0000 0.0000 0.0000 0.0000 1.3898 ' 

Dci = 0.0000 0.0000 ‘ 

- • J 

• Control augmentation compensator : Forward loop controller 

yji = Cfi 

Cf^ = 0.2 in the longitudinal axis and for the lateral ca.se it can not be found 
out as the scale given in the MIL specification is not sufficient. 

Flight condition : Forward speed of 200kph 


t Stability augmentation closed loop control law : 

As in the above case of hover q and 9 are used in feedback. The controller is 

■ \ B ^ 

Oi — : 

^ Cci -Dci ) 

where 


- 0.2829 

- 0.3710 

- 26.5326 

14.8712 

446..5495 

0.0891 

- 0.4112 

5.1943 

- 1 . 64.59 

- 77.0964 

= 0.0000 

0.6779 

- 27.7090 

12.8070 

432.6279 

0.0000 

0.0000 

28.0092 

- 14.7179 

- 458.3225 

0.0000 

u 

0.0000 

0.0000 

0.1179 

- 1.5189 

2.0253 

3.7937 




0.9013 

0.5828 




= 0.2947 

0.8890 




- 2.8761 

- 3.1318 




0.0984 

0.1166 





‘^7 


Cc, = ; 0.0000 0.0000 O.OOOO O.OOOO I.TGOo ! 

J 

= 0.0000 0.0000 I 

1 In* lal( lal (lyii«iini(s is slabh*. So ilicrr is no lunnl cjf dcsi^^ni I, In* <'on{.r<)tl(‘r. 

• Control augnientalioii compensator :Forward loop coiit roller 

C/j = 0.2 for the longitudinal model and again for the lateral model the scale 
is insufficient. 

3. Robust stabilization : 


• Nominal condition : Hover 

Multiplicative uncertainly M{Po,\V) is used where 


' (10s+l)(0.l3+lj 

( 152 + 2 . 84 + 1 )' 

: 0 


M{Po,\V) includes the full forward speed. 


0 : 

(0.25+2) ■ 

{2j+1)’ j 

Controller docs not e.xist in llii.s 


case. 


• Nominal condition : Hover 


ir 

V'c 

C 


(0.35+0.1) 

(0.15+1) 


C«c 


.4 B 


C D 


■T 

0 ‘ 

(0..35+0.1) 
(0.15+1) j 


.4, B, C are 

-40.2498 
-8.2440 
-0.4122 
-4 = 861.9700 

0.0000 
0.0000 
0.0000 


0.0055 0.0385 

-0.01577 0.0312 

0.0048 -0.3346 

0.9390 -0.0905 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 


42.8005 7.0336 

-3.2054 -3.2860 
36.9128 48.7944 

0.3285 0.3571 

1.4783 0.77-50 

0.3725 0.5790 

-0.0195 1.0569 


358.0952 

-(>7..S3!]i| 

•30.9651 

29.9135^ 

-356.2358 

— K)9.37(ij 
-3.4087 1 

-1.5.8961 

-4.7164 

1 

-3.9325 

-5.5(i42 ! 

0.1889 

-10.5!7m 



-0.1360 -0.0960 ■ 

0.0074 0.0181 

-0.0-563 -0.2.5-54 
-O.OOGl -0.004-5 
-0.0004 -0.0038 
-0.0002 -0.0034 
0.0000 - 0.0002 

-0.24980.00-55 0.0385 5.6181 -0.0884 0.0198 0.0792 
0.0000 0.0000 - 
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Chapter 4 

Resuhs Anc Discussion 


Seveial simulation studies were performed to analyse the dynamic performance of the 
closed loop system for controllers designed in the last chapter. The results of these studie.s 
are presented in three sections. The first and second sections describe the simulal.ioii 
results of the decoupled control laws developed at hover and at a forward speed of 200k])h. 
In the last section, the results pertaining to the simulation of the closed loop full order 
system at lio\'cr is presented. 

4.1 Simulation results of the decoupled control laws at 
hover 

The mathematical model of the helicopter at hover is given in Appendix D. The system 
dynamics is decoupled into the longitudinal and lateral axis components and two separate 
controllers are designed for each axis as described in section 3.4. 

Longitudinal dynamics 

The longitudinal dynamic model of the helicopter at hover is unstable and the control law 
of Eqn.(3.10) makes the closed loop system stable. The poles of the uncompensated and 
compensated longitudinal model is given in Table 4.1. From Table 4.1, it can be seen 
that the uncomjjensated longitudinal model has one unstable complex conjugate' pole. 
This pole can be relocated in the left half s-plane using pole placement techniques. In 
the compensated the poles have been relocated at -0.2±0.333ii. The sensitivity function 
S, of the compensated system is found out and |j 5 ||co has been calculated as 9.oC40. 
The high ^'alue of || 5' lU. indicates that the system is highly sensitive. The reason for 



Table 4.1: Pole.s-lougitiulinal clynaniics 


uncompensated system 
-1.4050 

0.0356 ± 0.3331/ 

-0.3386 

-40.0000 


compensated system I 
-1.4050(2) 

-0.2 ±0.3:331/(2) 
-0.3386(2) 

-40.0000(2) ; 


this high sensitivity may Ije due to the fact that sensitivity minimization has not been 
considered in the design criteria. The response of the closed loop system for unit step 
input in longitudinal cyclic(0ij) are shown in Fig.(4.1)-(4.4) for all the state variables. 
For comparison, the step response is evaluated for three different pole locations. From 
Fig. (4.1)' it can be seen that for a longitudinal input, the forward velocity [uh) of the 
vehicle will change to a new value. However it must be noted that the response curves 
are not valid beyond initial few(2-3) seconds. Because, once the vehicle starts moving 
forward the system transfer function will change. 

The step response of the vertical velocity component Wk is shown in Fig.(4.2). The 
change in lOh is very small. So it can be said that the influence of Ou on to/, is minimal. 
The response of pitch rate q to longitudinal cyclic input 0i, of uint magnitude is shown in 
Fig.(4.3). The response of pitch angle 0 to a unit step longitudinal cyclic input is shown 
in Fig.(4.4). Comparing the figures, it is observed that for the complex conjugate pole at 
-0.2 ± 0.3331/, the initial response is more sensitive in comparison to the other two pole 

locations(— 0.4 ± 0.3331/ and —0.7 ± 0.3331/). 


Lateral dynamics 

The lateral dynamics of the helicopter at hover is unstable. So a controller has been 
designed to make it stable. The poles of the uncompensated and compensated system 
are shown in Table 4.2. The response of the compensated system for a unit lateral cyclic 
step(«,.) are shown in Figs.(4.5H4.8). Each of the figures shows the respome of all the 
states, for three different pole locatioiis(- 0.2 ±0.33311. -0.4 ± 0..i.i.!lt. -O.i ±0..i331/). 
The results shown in Figs.(4..5H4.8) indicates that for the pole a, -0.2 ± 0.333,, is „,o,-e 
sensitive than at otlier locations. 



lable 4.2; Polc.s-laleral clynaiiiics 


niiconij)('nsa(('(l .systx'in coinix'iisaltnl svslt'in 

-5-21 SS -5.2188(2) ^ 

-0.8794 -2.0000 (2) ‘ 

0.0021 ± 0.3407? -0.2 ± 0.3.33h'(2) 

-40.0000 -40.0000(2) 

4.2 Simulation results of the decoupled control laws at 
200kpli 

A controller for stability is designed for helicopter operating at a forward speed of 200kpli. 
The poles of the uncompensated and compensated longitudinal dynamics is given in Tal)l<‘ 
4.3. The response of the compensated system to a unit step input in longitudinal cyclic 

Table 4.3: Poles- longitudinal dynamics at 200kph 

uncompensated system ' compensated system 
-40.0000 -40.0000(2) 

-0.4595 ± 0.3579?' -1.4050(2) 

0.6359 ± 0.4262?' -0.3386(2) : 

■ ^ -0.2±0.3331?(2) 

(^ll 5 ) is shown in Fig.(4.9) for all the variables. The characteristics of the initial response is 
similar to that observed in hover. The lateral d^mamicsof the model at 200kph is found to 
be stable. The poles of the lateral model are —0.6894 i 2.3547?, —5.8189, —0.0416, —40. 
So there is no need of designing the controller for stability augmentation. 





4.3 Simulation of the controller for the coupled systeiii 
at hover 

The coupled system behaviour is investigated using the decoupled longitudinal and lateral 
controllers. The poles of the coupled system both uncompensated and compensated are 
given in Table 4.4. 


Table 4.4: Poles-coupled system 


uncompensated system 
-40(2) 

-5.1082 

-1.3138 

-0.7953 ± 1.2689f 
0.1375 ±0.3462?: 
-0.0143 ±0.3609?’ 


compensated system 
-39.9937(2) 

-0.3348 

1.0781 

-40.0039 ± 0.0097?’ 
0.1744 ± 0.5198? 
-5.1944 ±0.1040?’ 

, -2.2096 ± 0.4830? 
-1.1687 ± 1. 549 1?’ 
-0.3405 ± 0.6973?’ 
-0.33410.4429? 
-0.1827 ±0.1150?: 


It can be seen from Table 4.4 that the uncompensated system is unstable. It does not 
satisfy MIL-H-8501A stability requirements because the frequency of oscillation of the 
unstable mode is 18.14 seconds. Its real part is such that the magnitude of the signal 
doubles in every 5.04 seconds. The compensated system with decoupled longitudinal 
and lateral controllers(designed for pole placement at -0.2 ± 0.3331?) is also unstable. 
In addition it does not meet the MIL specification on stability requirements, since the 
frequency of oscillation of unstable mode is 12.57 seconds with the real part such that its 
magnitude doubles in every 4.23 seconds. 

The response of the .system was calculated for two different controllers designed for 
placement of poles at either -0.2 ± 0.3:331i or -0.4 ± 0.3331i. It can be seen fron. the 


figures that for bot h set of controllers the couj)led system shows unstable behaviour. Tlie 
system with the controller designed for poles at — 0.4±0.3331i seems to be more unstabh' 
compared to the controller designed for poles at -0.2 ± 0.3331b 

Since the compensated system with the decoupled controllers is unstable, a new con- 
troller is designed for the full order system. The poles of the uncompensated and com- 
pensated system are shown in Table 4..5. Thus it can be seen from the Table that the 
compensated system is stable. 


Table 4.5: Poles-coupled system 


-5.1082 
-1.31.38 
0.79.53 ± 1.26891 
0.0143 ± 0.36091 
0.1375 ±0.34621 
-40.0000(2) 


- 5.0000 . 
- 2.0000 , 
-0.5 ± 1.26891 : 
-0.3 ± 0.41(2) 
-0.2 ±0.34621(2) ' 
-40.0000(4) 


The step response of the compensated system are shown in Fig.(4.1S) and (4.19) for lon- 
gitudinal and lateral state variables. From these figures it can l:)e seen that when both 
longitudinal and lateral cyclic step are appied to the compensated full order model, the 
longitudinal and lateral variables behave in a way similar to the decoupled longitudinal 
and lateral model. 

Based on the limited results presented above, one can conclude that care must be 
exercised in using the controller designed for the decoupled dynamics for the full order 
system. Hence it may be beneficial to design for the full order system. 
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Figure 4.1: Response of Uh to a step i/p in ^i^-Hover 
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Figure 4.2: Response of Wh to a step i/p in ^u-Hover 
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Figure 4.3; Response of q to a step i/p in ^i,-Hover 
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Figure 4.4: Response of 0 to a step i/p in ^?l 5 -Hovcl 
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Figure 4.5: Response of Vh to a unit step in ^le-Hover 
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Figure 4.6: Response of p to a unit step in 0ic-Hover 
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Figure 4.7; Response of r to a unit step in ^ic-Hover 
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Figure 4.8: Response of 0 to a unit step in ^ic-Hover 
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Figure 4.9: Response of all longitudinal state variables to a unit step in di,-200kph 
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Figure 4.10: Response of u/, to unit step in both and ^jc-Hover 
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I'igure -11 J; Response of Vh to unit step in both 9^ and ^ic-Hover 
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Figure 4.12: Response of lOh to unit step in both du and 0ic-Hover 
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Figure 4.13: Response of p to unit step in both 9 .nd 0ic-Hover 
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Figure 4.15: Response of r to unit step in both and 6'ic-Hover 
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Figure 4.16: Response of 6 to unit step in both and f?ic-Hover 
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Figure 4.17; Response of 9 to unit step in both and ^ic-Hover 
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Figure 4.18: Response of the longitudinal state variables to unit step in both and 
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Figure 4.19: Response of the lateral state variables to unit step in both and <?ic-IIover 



C lapTer 5 

Conc-ucing RemarL<s 

In this thesis, \v(> discnssccl the application of optimization techniques and factoriza- 
tion theory' to Uk' (h'sign of control laws for flight control system of a helicopter. In the 
first stage the characterization of all closed loop controllers which assigns the closed loo]) 
poles in tlu' region for good handling qualities was carried out. This was developed using 
the standard cont rolh'r i)arametrization of the factorization theory. All these controllers 
augment th«' stal)ilit y of the vehicle as per MlL-H-SoOlA specification. In the next step, 
the forward coiupcuisat-or to augment the pilot input response is designed. Using the 
design rnetliodology outlined in this chapter, two separate two degrees of freedom con- 
trollers arc designed to gi\'e <lccoupled longitudinal and lateral control of the helicopter 
for two different, flight conditions. The simulation study presented has shown that the 
control strategy ]n'oi)osed gives very good performance in the decoupled longitudinal and 
lateral axes. But. the wliole coupled system does not perform well with the decoupled 
longitudinal and lateral controller, s. Hence using the same design strategy closed loop 
controller i,s designc'd for the whole coupled .system. 

Robust closed loop st abilization of the wliole flight envelope, by including two extreme 
flight condit ions imhox'er and full forward speed is investigated and is shown that for the 
uucertaity model chosen such a controller did not exist. Robust stabilizing contiollci is 
next <le,signcd at. hoven- by considering the multiplicative description of plant uncertainty 
which is more realistic. I'he state space algorithm of [8] is used in designing the robust 
stabilizing controller. 


5.1 Scope of future work 


• One of tlie short coinings of the designs so far undertaken has been that they have 
been restricted to one operating point. It is intended to develop a control system 
to cover wider operating envelope and this will entail undertaking designs at a 
number of operating points and investigating ways in which they may be scheduled 
or switched in a bumpless fashion. 

• A very important area is that of actuator saturation, which can lead to a loss in the 
direction of control input signals and thus to a degradation in performance. In thi.s 
tliesis we have not considered the problem of actuator saturation. It would be very 
important to develop a controller conditioning strategy to ensure that performance 
is maintained in the face of actuator saturation. 

• In this thesis, the controllers obtained are of high order(ie, 5 and above). The 
modern hardware does not pose serious limitations to implementing high order 
controllers. How(>ver it is still important to find reduced order controllers and 
confirming the performance of the closed loop system after reduction, due to their 
ease of implementation. 



APPENDIX A 


Parametrization Of All Stabilising Controllers 

In this appcnclix the formula for obtaining a doubly coprinie stable factorization of 
a plant transfer matrix P, is presented. Using this factorization, the formula for all 
stabilising compensators which guarantee poles of the closed loop system in a specified 
region of the complex plane, is also presented. 

Consider the system transfer matrix P represented in a state space model as 

.r = .4.r + Bu 
y = Cx+Du 

which is denoted in shorthand as 

^ A B 

P = . =C{sl -A)-^B + D 

The doubly coprime factorization of P is defined as follows. 

Definition:- 

A doubly coprime factorization(DCF) of P over a closed region fi in the complex plane 
is a factorization of P as 

P = iYp-i = 

such that the following properties are satisfied. 

1. N, D, Ah D have no poles in U. 

2. There are transfer matrices X, Y, .Y. V which have no poles in 0. 

3. Y, V' are square and nonsingular.' 

4. N, D, -A'. D. A', Y, .Y. Y satisfy the identity 

y A D -X _ ^ 

-iV D . N Y 
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The procedure for obtaining N, D, N, D of a given plant P from its state space discription 
is given below. Suppose the given system 

( .1 n ^ 

p = 

yc DJ 

is minimal. Select matrices K and H such that Ak •'= A — BK and Ah := — HC have 

no eigen values in Cl. Then N. D, Ah D, X. Y\ A", Y are given as [14] 

/ Ah B- Hd'' 
yC D ! 

' .Ah H ^ 

\-C I ! 

( .-Ik B '' 

C - DK D I 

f .Ik B '' 

{-h- I J 

f .Ah iP 

\ 0 1 

( .Ah B-Hd\ 

\K I , 

Pak iP 

V 0 J 

< ..Ik II '' 

yC-DK I j 

Formula For All Stabilising Controllers 

In the DCF above, if Cl is a compliment of an open connected domain P in the left half 
of the complex plane, then the factors N, D, Ah D, A', 1’, .A, Y are all stable traush'r 
functions with all the poles in F. The well known formula which makes use of the DCf 
and expresses all the stabilising compensators C(s) such that the closed loop system of 
Fig. (5.1) with the transfer matrix from w to z has all its poles in F is given as follows. 
Let Q be any proper rational transfer matrix with no poles in H then 

C = {[Y-QN){X + QDrK det{Y-QN)^0] 

f\7 


N = 

D = 

A' = 

D = 

X = 

T = 

A' = 

Y = 



Figure 5.1: The Standard Block Diagram 


= { {X + DQ)iY - XQr\ det{Y -NQ)^0} 

It is well known from the stable coprime factorization theory [14], that if C is a stabil- 
ising controller expressed as above then the closed loop transfer matrix of the system in 
Fig.(5.1) can be expressed in terms of the DCF and Q, where Q appears in the transfer 
matrix in an affine manner. This mathematical fact turned out to be of vital importance 
in the development of modern feedback design theory. 


fi8 



APPENDIX B 


The Four Block Control Problem 

'I'lic ,s(aii(lar<l l-hlock problem is shown in the figure below. In the modern four bitjck 
ai)proaeli to control system design the plant G, to be controlled is assumed to have the 
structure as shown in Fig. (5. 2). 

G = 

2 “ > *. 


Figure 5.2: The Generalized Plant 

where w is the e.xogeneous in])ut, typically consisting of command signals, disturbance's 
and sensor noises, u is the control signal, z is the output to be controlled, its components 
typically being tracking errors, filtered actuator signals etc. y is the measured output. 
A feedback controller (.’ is attached to the plant to form the feedback system shown in 

% 

T 


Figure 5.3: The Standard Problem 


Fig.(5.3). 


vJ 


~ > 






For G given l)y tlu' four block Iransfcr matrix 

r 

^MI G'i2 

. ^ Gbi G-22 

The equations of (he feedback system are then given by 

c = G’liU' + G'nn 

y = G^iJc + Gbau 

u — —Cy 

Hence for a well posed system(ie, when det(I+GC)^ 0), the transfer matrix denoting the* 
response z due to exogenous input w is obtained as 

y = 

where 

r^u, = Gii — GnC{l + G22G) ^Gbi 

Standar'd problem 

The standard problem is to find a real rational proper C to minimize the //oo norm 
of T:,!,, the transfer matrix from w to z under the consraint that C stabilizes G. All 
the requirements like robustness, disturbance rejection, tracking accuracy etc can be 
captured in the transfer matrix between w and z. The standard problem above is quite 
well developed and a state space algorithmic procedure is available for its solution. Oik' 
such procedure is detailed below. 

Glover-Doyle state space algorithm 
The 4-block plant G is considered in the state space form 

^ A' Bx B 2 ^ 

G{s) = Gi Dll Di2 ' 

^ G'2 .D21 D22 j 

with the following assuptions. 

1. [AJh) if^ .stabilizable. (G2,.4) is detectable. 
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2. D\i is full rolmnn rnnk, Du is full row rank, 
.4 - ju.-/ B2 


•3. 


Cl D 


has full coloumn rank Vo;. 


12 


.4 - io.’/ Bi 

4. has full row rank Vtu. 

C, D21 _ 

Define the two Hamiltonian matrices as 

A 0 


: 

•^^(7) = 


B 


-CfCi -AT' _ ^-CfDi. 

A'T' 0 CT' 

-BiBf -A ' -BiDT 


BR-^ DjCi B^ 


R~^ D.iBf C 


where 


Di. = 
'd.i = 


On On 
On ' 

Dll 


and 


R = DTDu- 


T 


R = D.iD\ 


I'lml 0 

0 0 

J 

tU,, 0 ' 
0 0 


Dn is partitioned as 


D 


11 


Dun Dim 
Dll'll D1122 

Let A, Q, R. be real nxn matrices with Q and R s\’mmetric. Define the 2iix2n Hamiltonian 
matrix 


H 


A R 
Q - 4 ^_ 

The dom(R.ic) consists of Hamiltonian matric(\s witli two proiierlic’s, stability ])ropert.y 
and complementarity property [13]. The conditions for H to belong to doni(Ric) are as 


r I 



follows. H should have no imaginary eigen values. R should be either positive semidefinit.e 
or negative semidelinite and (A,R) stabilizable. 

Theorni:- 

There exists stabilizing controllers such that || 7Ti. ||< 7 if and only if the following 
conditions hold: 

1. 7 > 777 «.r (cr[L>nn Dnvi]^ ^[D]n\ DJui]) 

2. e dom(Ric) and X{j) :=Ric[//.^] > 0 
Jocil') € dom(Ric) and >''(7) :=Ric[.y.x.] > 0 

3. p[A'(7)>-(7)] < 

When these conditions hold, one such controller is 

R-i 

A'(U = , . 

Ve Dj 

where 

b = -Dini Dn22 

C = [F2-I)(C2 + Fi2)]Z 
B = -H2 + {B2 + Hi2)b 

.-i = A + HC + [B 2 + Hn)C 

Z := (/-7-’^y'A)-^ 

p'' = [^/i 4 = -(A'/y + c’/'/;n)/r» 

H = [Hn Hi 2 H 2 ] = -(VC''’ + 



APPENDIX C 


Robust Stabilization 


The robu.t, stability condition, in teens of the standard 4-block |>t„ble,n te<,nin. 
II ||»< 1. The solvability of this standard problen. using ...odel rnatcl.mg tl.eory , 


given below. 

The transfer matrix between w and z can be written as [6] 


r^u, = 1\ - TiQT:,.. 


]• \ T Ts -iiul I\ arc defined as follows. Let the slandai 
with C(s) given in .Appendix -A. /i- ^2 And i.i aic 

set up of Fig.(5.1) in A,)|>e.Klis .A has stale s|.ace realisation 

( .4 Ih Ih ^ 

G{s) = Cl Dll Di 2 

C 2 D'li D 22 j 

Choose F and H such that dr := d - BF and d„ := d'- UC are stable. Then V„ 


and 2s are given by 


2’,(s) 


^ A R ^ 
{ C Dn ) 


Af IhF 

zi = 

0 Ah 

^ a 

Bi 

B = 

Bi-llD 2 i_ 

a = Ci-DnF D12F _ 

( Af B 2 ^ 

Cl - D12F Di 2 j 


Us) 


^ Ah Bi~BD2i^ 
('2 D21 j 



L(’l i,, / “■ 1. 2, (leiiolt’s the right lialf phiiie tran.snil.s.sloii zercje.s of 72 and 

Tlien inorder for tlie roljust stabilization ]>roblem to be solvable, Ti must satisfy IJk* 
condition, || T\{~i) || < 1 V? = 1,- 2, 
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APPENDIX D 


Model of The Helicopter at Different Flight Conditions 
The model of the helicopter at hover, 200kph and 330kph are given below. 


Model at Ho\’er 

System Matri.x (8X8) 


-0.1567£- 01 

0.1900E- 02 

0.3119E- 01 

-0.7445E - 02 

0.5274E - 02 

-0.1587£-02 

-0.4785E - 1 1 

-0.1703£ + 00 

1 

i 

1 

-0.5S33£' - 02 

~0.0007£ - 01 

0.1331E+0I 

-0.8362£ - 02 

0.6488E-0I 

0.9020E - 02 

0.170I£ + 00 

0.8476E - 03 ; 

, n 

7 

0.4771E- 02 

-0.5832E - 02 

-0.3346E + 00 

-0.1876E - 03 

0.4477E - 02 

0.1079E-01 

0.8369E - 02 

-O.I723£ -. oll| 

1 

•«0.2627£+ 01 
A = 

0.9390E + 00 

~0.0690£+01 

-0.30I4E + 02 

-0.5I62£ + 0i 

-0.253l£-f 01 

-0.426I£-0I 

-03097£ - 07 

-0.21 15£- 00 


-0.478 1£ + 00 

-0.9053E - 01 

0.1959E + 00 

-0.1322E + 01 

0.532DE - 01 

~0.9409£ - 09 

-0.1 HOE - 07 


-0.2003E + 00 

0.5671E + 01 

-O.I I70£ + 03 

-0.5650E + 00 

-O.6537E+0I 

-0.93I9£ + 00 

-0.i935£- 07 

-0.53ii£-;^ 

1 

O.OOOOE + 00 

O.OOOOE -f 00 

O.OOOOE + 00 

O.IOOOE+Ol 

-0.4977E - 02 

OJOnE-f 00 

O.OOOOE 4- 00 

O.OOOOE + m 


^ O.OOOOE + 00 

O.OOOOE + 00 

O.OOOOE + 00 

■ O.OOOOE + 00 

0.9900E + 00 

0.4915E- 01 

O.OOOOE + 00 

O.OOOOE + 00 

' 


Control Matri.x 

0.2264264£ - 0 1 -0.2060993S' + 00 

0.1950103£' + 00 0.1500457£ - 01 

-0.3645610£' - 01 -0.1030504S: - 01 
0.0306200£' + 02 -0.2353 194£' + 01 

B = 

0.1333310S + 01 0.2154925£: + 02 

0.1430312£’ + 02 -0.372050 !£ + 00 

0.0000000£ + 00 O.OOOOOOOS + 00 
O.OOOOOOOE + 00 O.OOOOOOOS' + 00 


Model at 200kph 

System Matrix (8X8) 


-0.4132E- 01 

0.5451- 03 

0,2400£-01 

-0.0I20E- 02 

0.3598E-01 

-O.96Q0E - 04 

-0.225 1£- 11 

-0.17I1E + 00 

-0.7167E- 03 

-0.2313E + 00 

-0.2185E- 01 

-0.4072E - 01 

-0.5673. £ - 02 

-0.9036£ + 00 

0.1716E + 00 

1 

O 

1 

0.6574E-01 

-0.i639£- 01 

-O.I360E + 01 

-0.2293£ - 01 

0.9500E + 00 

0.1633£- 01 

0.4800E - 02 

0.474 0£ - 02 

^ -0.1051£ + 01 

: 0.9650E + 00 

-0.3872E + 01 

0.1049E + 01 

-0.5592E + 01 

0.8396E + 00 

-03361 E - 00 

-0.1732E - 07 

-0.6403E - 09 

-0.5012E + 00 

-0.2310E-01 

0.I676E + 00 

Q.1754E+QI 

0.5943£ - 01 

-0.1865E- 00 

-0JU7E- 07 

-0.3683E + 00 

0.5887E + 01 

-0.3427E + 00 

-0.1066E + 01 

0.3539E- 01 

-0.I416E + 01 

-0.1025E- 07 

-0.5384E - 09 

O.OOOOE + 00 

O.OOOOE + 00 

O.OmOE + 00 

O.IOOOE + OI 

0.79 19£ - 03 

-0.2770E + 01 

0.0000£ 4* 00 

0.0(K)0£ 4 - 00 

0.0000£-h00 

O.OOOOE + 00 

O.OOOOE + 00 

O.OOOOE + 00 

Q.9996E + 00 

0,2857£- 01 

O.OOOOE 4 * 00 

0.0000£ 4- 00 


(0 


Control Matrix 


~-0.00.10756/r - 02 

--0.171 9079£ + 00 ' 


1 

0,I989993£ + 00 

-0.1195507£- 01 


l' 

' 0.19556I8E+ 00 

-0.1342436£ + 01 


f 

0.9I26846£' + 02 

B ~ 

0.162401 !£■ + 00 

-0.0747376£ + 00 

0.1937637£+ 02 

. 

h 

0.1954920£' + 02 

-0.1055517£ + 01 ; 


I 

0.0000000£ + 00 

0.0000000£ + 00 , 


1 

0.0000000£ + 00 

0.0000000£ + 00 




Model at 330kph 

System Matrix (8X8) 


-0.6001£-01 
0.4894£ - 02 
0.1189ir + 00 
-0.9696£' + 00 
0.5496i: + 00 
-0.1450£: + 01 
O.OOOO^ + 00 
O.OOOOE + 00 


0.705 1£ - 02 
-0.3597£ + 00 
-0.5054E-0I 
-0.731 1£ + 01 
-0.I562£ + 00 
0.7075£ + 01 
0 . 0000 £ + 00 
0 . 0000 £ + 00 


O.OlOOE-01 
-0,1541£ -01 
-0.14^9£ + 01 
0.5599£ + 01 
0,5074£ + 0I 
0.3160£ + 01 
0 . 0000 £ + 00 
O.OOOOE + 00 


-0.4801£ + 02 
-0.1150£ +00 
-0.4705£ - 01 
-0.4151£ +01 
0.5549£ + 00 
-0.4575£ + 00 
0 . 1000 £ + 01 
O.OOOOE + 00 


0.11I0£ +00 
-0.0503£ - 02 
0.1583£ + 01 
-0.1150£ + 01 
0.2164£ + 01 
0.4215£ -01 
0.6211£-02 
0.9903£ + 00 


0.3110£ - 02 
-0.1577£ + 01 
0.2015£-07 
-0.1502£ + 00 
0.7040£ - 01 
-0.2052£ + 01 
0.I063£ + 00 
0.5034£-01 


-0.4326£ - 1 1 
0.1699£ + 00 
0.9930£ - 02 
-0.201 1£- 07 
-0.1543£- 06 
-0.1496£ - 07 
0 . 0000 £ + 00 
0 . 0000 £ + 00 


-0.1702£ + 00 
-0.1057£ - 02 
0.1009£- 01 
-0.923 1 £ - 09 
-0.I098E- 07 
0.1255£ - 09 
0.0CK)0£ + 00 
0 . 0000 £ + 00 


Control Matrix 


-0.4634399£ - 02 
0.2312920£ + 00 
0.418931 1£ + 00 
0.838 1 099 £ + 02 
-0.1757970£ + 01 
0.1294223£ + 02 
O.OOOOO0O£ + 00 
0.0000000£ + 00 


-0.1160004£ + 00 
0.1221094£-01 
-0.2252474£ + 01 
0.1546433£ + 02 
0.3070232£ + 02 
0.47475 12£ + 01 
0.0000000£ + 00 ' 
0.0000000£ + 00 
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